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Differences between lattice Monte Carlo data and perturbation theory (for example the lack of asymptotic 
scaling) are usually associated with the 'bad' behaviour of the bare lattice coupling go due to the effects of large 
(and unknown) higher order terms in go- In this philosophy a new, renormalised coupling g is defined with the 
aim of making the higher order coefficients of the perturbative series in g' as small as possible. 

In this paper an alternative scenario is discussed where lattice artifacts are proposed as the cause of the 
disagreement between Monte Carlo data and the go-perturbative series. We find that with the addition of a 
lattice artifact term, the usual asymptotic scaling expression in go is in excellent agreement with Monte Carlo 
data. Lattice data studied includes the string tension, the hadronic scale ro, the discrete beta function, M p , f n 
and the IP-IS splitting in charmonium. 



1. Introduction 



2. Lattice Distorted-Perturbation Theory 



A necessary condition for lattice predictions 
of QCD and other asymptotically free theories 
to have physical (continuum) relevance is that 
they reproduce weak coupling perturbation the- 
ory (PT) in the limit of the bare coupling go — > 
0. This perturbative scaling behaviour (a.k.a. 
asymptotic scaling) has not yet been observed for 
complicated theories like QCD when the bare lat- 
tice coupling go is used as the expansion param- 
eter. 

As a result of this disappointing disagreement, 
various workers have proposed methods of im- 
proving the convergence of the perturbation series 
by a re-expansion in terms of some new coupling 



This paper studies an alternative viewpoint in 
which the disagreement stems from lattice arti- 
facts f|. In this talk, it is shown that these 
terms can provide the mismatch between the lat- 
tice Monte Carlo data and go~PT without resort- 
ing to the use of a re-defined coupling g' . 

The QCD quantities studied in this analysis 
are: the string tension, ^/a; the hadronic scale, rn 
0; M p \ fn', the IP — IS* splitting in charmonium; 
and the discrete beta function A/3. 

The results discussed here are presented in 
greater detail in [||. 



Two-loop perturbation theory predicts the run- 
ning of the lattice spacing a with coupling g 2 as 
follows, 
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fpT(g 2 Y 



where 
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Lattice calculations predict a by calculating 
some dimensionful quantity on the lattice, and 
comparing it with it's experimental value. As 
is well known these lattice values do not follow 
the above perturbative behaviour (when the bare 
coupling go is used). There are a number of pos- 
sible causes of the disagreement: quenching; fi- 
nite volume effects; unphysically large value of the 
quark mass; a real non-perturbative effect; the in- 
clusion of only a finite number or terms (i.e. two) 
in the PT expansion; and lattice artifacts due to 
the fmitcness of a. For the reasons outlined in 
||, the first three effects cannot give rise to the 
sizeable discrepancy between lattice data and PT. 
As far as true (i.e. continuum) non-perturbative 
effects are concerned, the overwhelming expecta- 
tion is that for cut-offs of a -1 <; 2 GeV these ef- 
fects should be minimal. Therefore the disagree- 
ment can only be due to either or both of the last 
two possibilities. 



2 



Fitting 








a ^ from 








Method 








M p 


u 


IP - IS 


A/303) 


<7o-PT 


A [MeV] 
X 2 /dof 


1.254(3) 
484 


1.599(5) 
262 


1.63(1) 
10 


1.48(2) 
9 


1.45(5) 
6 


702 


Leading- 
Order 
LDPT 


\ r at 

A [MeV] 
X 2 /dof 


0.204(2) 

1.90(1) 

3 


0.150(2) 
1.958(9) 
16 


0.22(2) 
2.15(5) 
1.1 


0.34(3) 

2.2(1) 

1.6 


0.35(6) 

2.5(3) 

0.3 


0.373(5) 
4.2 


Next-to- 
Leading- 
Order 
LDPT 


VOL 

A n 

VOL 

A n+2 

A [MeV] 
X 2 /dof 


0.26(2) 
-0.024(6) 
1.96(2) 
1.7 


0.29(1) 
-0.046(3) 
2.14(2) 
1.4 








0.24(1) 
0.050(5) 

1.7 




A [MeV] 
X 2 /dof 


17.34(4) 
160 


20.89(7) 
47 


21.4(2) 
1.3 


21.0(3) 
2.5 


19.3(7) 
1.5 


78 


9E-PT 


A [MeV] 

*7*>/ 


4.81(1) 
52 


5.56(2) 
15 


5.77(4) 
3.6 


5.58(7) 
1.4 


5.2(2) 
0.3 


19 



In this section we study the effect of lattice ar- 
tifacts. These can be parametrised (to leading 
order) by modifying eq.(P as follows: 
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(with no implicit summation over a or n). 

Note that the 0(a n ) coefficient in Eq(g) has 
been normalised so that X" is the fractional 
amount of the 0(jjq a n ) correction at a standard 
value of go = 1 (i.e./? = 6). For M p ,f v and the 
IP — IS splitting, a — & n = 1; for a and A/3, 
a = & n = 2; and for r , a = 2 & n = 2. 

Lattice Monte Carlo data taken from many dif- 
ferent collaborations are fit to eq.(||) (see || for a 
list of these references). This provides the values 
for A, X™, and the \ 2 as listed in the table. Also 
shown in the table are the fits to (2-loop) go - PT. 
(The fit in this case is eq.(0) with g = g .) 

We see that leading order "lattice distorted- 
PT" fits the data very well compared to the go- 
PT case, with the \ 2 /dof down by an order of 
magnitude or more. 

As a further check of the method, we include 
in the fit the next-to-leading term in a. How- 
ever, due to the large statistical errors in some 
of the lattice data we perform this fit only for 
cr, ro and A/? where the statistical errors are very 
small. The fitting function appropriate for these 



quantities is: 
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The results of these fits are displayed in the 
table. 

Obviously in the limit of infinite statistical pre- 
cision, adequate fits to the lattice distorted-PT 
formula would only be obtained if the 0(a n ) 
terms were included to all orders. The fact that it 
is necessary to go to next-to-leading order for the 
cr, ro and A/3 data to obtain a sensible x 2 /dof 
simply states that these data have sufficiently 
small statistical errors to warrant this order fit. 

The rest of this section comments on the results 
of these fits. 

It is clear that for a, ro and A/3 data, the agree- 
ment between the data and lattice distorted-PT is 
remarkable considering the tiny statistical errors 
in the lattice data. 

Another important finding is that the values of 
A for the various quantities are all consistent with 
A = 2.15 MeV within around la with the only 
exception being the string tension. This slight 
discrepancy can easily be explained by the ef- 
fects of quenching, and the uncertainties in the 
experimental value of a. Taking A = 2.15 ± 10% 
MeV as an overall average, and converting to the 
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MS scheme, we have A^£" 0) = 190 ± 20MeV. 

.MS 

This compares well with other lattice determina- 
tions and therefore supports the validity of this 
approach. 

The typical values of X™ in the table are 20- 
40%. A study at /3 = 6.4 || found that non- 
perturbative determinations of the renormalisa- 
tion constant of the local vector current vary by 
10-20% depending on the matrix element used. 
Since the spread in Zy en has been interpreted as 
0{a) effects fit], we can assume that O(a) effects 
of around 20-40% are reasonable at ft = 6.0. 

The coefficients for the second order terms are 
an order of magnitude less than the first order 
terms. This follows our expectation that the ex- 
pansion in Jpt in cq. (j^) forms a convergent series. 

One of the most exciting features of the lat- 
tice distorted-PT approach is that it can repro- 
duce the behaviour of A/3. The interpretation of 
the well-known discrepancy between g^-VT and 
Monte Carlo A/3 data has been problematic in the 
past. For example, in ||, a fit was attempted to 
their A/3 data using a coupling with two free pa- 
rameters. A good fit was obtained only for an un- 
physical value of one of these parameters, leaving 
the explanation of the discrepancy open. The lat- 
tice distorted-PT approach solves this problem. 

Finally, as far as the fit to M p , f n and the 
IP — IS splitting are concerned, the errors in 
the lattice data are large enough to allow many 
functional forms. Thus these data do not con- 
strain the lattice distorted-PT fit (or fits from 
other schemes). 

3. Fits Using a Renormalised Perturbation 
Theory 

In this section we fit the Monte Carlo data for 
a -1 to the following functional form: 
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where g' is some "renormalised" coupling which 
is in turn a function of the bare coupling go- Note 
that in this philosophy, the failure of asymptotic 
(i.e. perturbative) scaling is explained by higher 
order terms in perturbation theory and lattice ar- 
tifacts are assumed to be negligible. 



We studied two definitions of g 1 : (i) The MS- 
like coupling ||, 

1 
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, = \ < \TrU p i aq > MC +0.025. 



MS 

and the scheme of 



based on the plaquette, 
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1/3 



1- < 



^Trllplaq >MC 

(Note in the full paper, alternative definitions 
of g' are studied as well as those above IBJ.) 

The results of fits using these definitions of g' 
in the fitting function Eq. (^) are displayed in the 
table in the rows headed gj-fg and gE- 

As can be seen they do not reproduce the 
Monte Carlo results nearly as well as the fits from 
the lattice distorted-PT. 

4. Conclusions 

This talk studies the question of why (dimen- 
sionful) lattice Monte Carlo quantities do not fol- 
low the predictions of 2-loop perturbation theory 
in the bare coupling. The conventional answer 
to this problem is that higher order terms in g 2 
spoil the behaviour of perturbation theory, and 
that therefore an improved coupling is required. 
An alternative approach is presented here where 
the effects of 0{a) are shown to be able to provide 
the mismatch. The quality of the fits using this 
latter approach, and various arguments outlined 
in SecJ| support this philosophy. Further studies 
are required to unambiguously confirm this issue. 
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